load paramsDirectlyModified2dev.mat;
load paramsFixed.mat;

pStr2Double = defParamSet(paramsDirectlyModified2dev);

resParams0 = CalcResParams(paramsFixed, pStr2Double, 0);    
optsol0 = NNsysExact_dev(10000, resParams0);

resParams15 = CalcResParams(paramsFixed, pStr2Double, 1.5e+6);    
optsol15 = NNsysExact_dev(365, resParams0);

resParams75 = CalcResParams(paramsFixed, pStr2Double, 7.5e+6);    
optsol75 = NNsysExact_dev(365, resParams0);

save('otpsol.mat', 'optsol0', 'optsol15', 'optsol75');